Tunable thermal emission at infrared frequencies via tungsten gratings 
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The author investigates the manipulation of thermal emission by using one-dimensional tungsten 
gratings with different groove depths, ft is found that, by systematically increasing the depth of 
the groove, the linearly polarized emission at particular frequencies can be substantially enhanced 
to achieve that of the blackbody radiation limit, whereas the emission in other frequency ranges 
shows no noticeable changes. The results can provide useful insights into the design of thermovoltaic 
applications. 
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Thermal radiation heat transfer is a ubiquitous pro- 
cess in nature, and thermal emission has been an ac- 
tive research area for many years [l], 0. The main 
idea behind controUing the emission spectrum of an 
object is to manipulate its photon density of states 
(DOS). The photon DOS of a natural substance such 
as a metal or a dielectric, can be modified slightly 
by changing the temperature. The recent invention of 
man-made structures such as photonic crystals (PCs) 
and plasmonic surfaces (_PSs) has opened up the area 
of photon DOS control [1, i, i, @, [tI, i, S 13. Re- 
cent years have witnessed rapid growth in the study of 
thermal radiation control via such artificial structures 
[H [H, [11 B [S E [S [H, [S S ill. For exam- 
ple, radiation at infrared frequencies can be suppressed 
|14| by manipulating the geometry and material param- 
eters of metallic PCs with large band gaps over the in- 
frared frequencies, thus allowing radiation at optical fre- 
quencies to be more efficient. However, the complex- 
ity of the fabrication process for realistic PCs remains 
a challenge. A structure that can be easily fabricated 
while providing the fiexible control of thermal radiation 
is therefore highly desirable. Although the fabrication of 
bulk PC structures still presents a major challenge, tex- 
tured surfaces can be made relatively easily even for sub- 
wavelength patterns. In addition, the advantages of bulk 
PCs over surfaces with man-made structures in thermal 
radiation control are still being debated. For instance, 
it is difficult to fabricate the woodpile structure used 
by Fleming et al. and Han et al. [l^ Although the 
proposed use of self-assembly techniques [H f23| seems 
promising, the tungsten inverse opal structure [2J| fab- 
ricated by such a technique has not been shown to have 
substantial superiority over other known PC structures 
in controlling thermal emission. Although extensive at- 
tention is being focused on bulk PC structures, the mod- 
ification of thermal emission via textured surfaces has 
also been widely studied. Recently, Fleming et al. [13], 
Wan and Chan [l^ , and Yannopapas [1^ reported inde- 



pendently that a thin photonic crystal slab is sufficient 
to achieve the strong emission of electromagnetic radi- 
ation at desired frequencies, and that emission at other 
frequencies can be suppressed. Moreover, the thermal 
radiation effect of textured surfaces was also studied by 
Pendry Joulain et al. [H], and Fu et al. In 
view of these recent findings, textured surfaces appear to 
be promising candidates for thermal radiation control. 
In this study, the author uses a model system consisting 
of tungsten gratings and investigates the effects of groove 
depth on emission properties. It is found that, by system- 
atically increasing the height of the groove, the emission 
at particular frequencies can be tuned to achieve that 
of the theoretical limit, that is, the blackbody radiation 
limit, wheres the emission in other frequency ranges has 
no noticeable changes. Therefore, the tungsten gratings 
may be useful to the design of thermal photovoltaic de- 
vices, which require a structure that emits strongly at a 
particular frequency. 



h 
d 




an 



'Electronic address: |jwan@phy.cuhk.edu.hk| 



FIG. 1: Cross-session of the tungsten grating. The thickness 
(d) is equal to 0.5 /^m, and the other parameters are defined 
in the text. 

The cross-section of the model tungsten grating is 
shown in Fig.[T] The grating has a period of ao = 1.0 /im 
and a thickness of d — 0.5 /im. It is considered to be opti- 
cally thick because it is much thicker than the typical skin 
depth (~ 20 nm) of metals, thus allowing no light to pen- 
etrate it. The width of the groove is fixed at 0.3 /im, and 
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the depths are fixed at h = 0.5 and 2.0 ^m. According to 
KirchhofF's law, the spectral emittance ekC^) of a surface 
is equal to its spectral absorbance A^]f^{Lu) at thermal 
equilibrium; here, k and uj denote the wavevector and 
angular frequency of the emitted radiation. The validity 
of KirchhofF's law for textured surfaces has been demon- 
strated by Luo et al. [28|. Therefore, to understand 
the emission properties of tungsten gratings, the angle- 
dependent absorbances are calculated using the transfer 
matrix method (TMM) [2^. The absorbance for each k 
and oj is given by A_y^{oj) = 1 — i?k(^) — ?k(ci;), where 
R]i{uj) and T'k(w) are the total reflectance and transmit- 
tance of an incident wave. In our calculations, T^^io) = 
because of the large thickness of the gratings. 
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FIG. 2: Calculated (a) TE and (b) TM angle-dependent 
absorbance spectra for a grating of h = 0.5 jj,m. The high level 
of absorbance at E > 0.9 eV is mainly due to the intrinsic 
absorption of bulk tungsten. Surface plasmon excitation [Eq. 
([T])] in TE polarization is indicated by the green dashed lines. 
The dielectric data of tungsten is taken from Ref. 31. 
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FIG. 3: Calculated angle-dependent absorbance spectra for a 
grating of h — 2.0 /im. Note the absorption peaks ai E ^ 0.39 
and 0.62 eV, which are due to the excitation of the waveguide 
mode. 

The resultant absorbance spectra for TE (E _L y) and 
TM (H ± y) polarizations are shown in Figs. [2] and [S] 
For the latter, no noticeable difference in absorbance be- 
tween the different structures can be observed, and the 
overall absorbance of both structures is close to that of a 
plain tungsten slab. The strong absorption for E > 0.9 
eV (A ^ 1.38 /xmjis due to the intrinsic optical prop- 
erty of tungsten [30|. Therefore, the grating geometry 
does not significantly enhance absorbance for TM polar- 
ization. In contrast, the absorbance for TE polarization 
changes considerably with groove depth. In addition to 
the strong tungsten absorption at i? > 0.9 eV, strong 
absorption peaks are also observed for TE polarization 



at > 1.2 eV. The locations of the absorption peaks 
appear consistently and systematically for both h — 0.5 
/im and h = 2.0 /im. In addition, the peaks exhibit typ- 
ical band-folding behavior, thus indicating that they are 
due to the excitation of the propagating surface plasmon 
polaritons (SPPs) by the incident waves. To trace the 
origin of the absorption, the dispersion relation of the 
propagating SPPs is plotted on top of the absorbance 
spectra 
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where n is an integer, lo is the angular frequency, k is the 
in-plane wave vector in the a;-direction, which is related 
to the incident angle {9) by k = (uj/c) sin 6*, c is the speed 
of light, and Cair and ew are, respectively, the dielectric 
functions of air and tungsten [sT] . As can be seen, equa- 
tion (P) fits very well with the TE-absorption peaks for 
each grating, which indicates that the physical origin of 
the absorption peaks for E > 0.9 eV is the excitation of 
the propagating SPPs. 

We next focus on the region of E < 0.9 eV, in which 
absorption should be weak for tungsten. For h — 0.5 /im, 
non-dispersive absorption peaks are observed at ~ 0.45 
eV (A ^ 2.8 iim). As revealed later, the enhancement of 
absorption is by almost a factor of 10 compared to that 
of a tungsten slab. In addition, when the depth increases 
from 0.5 /im to 2.0 /im, multiple peaks are observed at 
E = 0.39 eV (A - 3.2 /im) and E = 0.62 eV (A ~ 
2.0 /im). Although not shown in this work, a similar 
scenario occurs if the groove width increases, whereas no 
considerable absorption peaks are observed for a groove 
width of less than 0.2 /im. Moreover, it should be noted 
that a,t E ^ 0.62 eV, the absorbance is very close to 1 
and remains strong even when the incident angle is close 
to 90°. These absorption peaks could be attributed to 
waveguide modes, as proposed by Porto et al. [s^, in 
which a standing wave is formed inside the groove and 
results in energy accumulation. 

To examine the dependence of TE absorbance on 
groove depth, the absorbance spectra for each depth at 
two incident angles are plotted (0° and 60°) in Fig. U) 
The absorbance of a thick tungsten slab is also drawn 
for comparison. As can be seen, the absorbance of a flat 
tungsten slab is low {A < 0.02) for near infrared (IR) 
wavelengths (A > 2 /im), which is due to the large dielec- 
tric mismatch between tungsten and air. As the mag- 
nitude of the real part of the dielectric function drops 
significantly for A < 1.4 /im, the dielectric mismatch is 
reduced, thus resulting in the observed increase of ab- 
sorbance for shorter wavelengths, and A reaches a plateau 
value of A ~ 0.6 for A < 0.8 /im. Therefore, the TM 
absorbances shown in Figs. [2] and [3] could mainly be at- 
tributed to the intrinsic absorption of tungsten. 

In the case of the TE mode, the effects of the groove 
can be clearly seen. Depending on different geometries 
and incident angles, absorbance can be enhanced by as 
much as a factor of 2 with respect to that of a tungsten 
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slab. For near IR wavelengths, the groove results in the 
great enhancement of absorbance at certain frequencies. 
For example, at normal incident {9 = 0°), A is enhanced 
by almost an order of 10 at A ~ 2.8 /j,m for h — 0.5 
fim, which corresponds to the peak at E ^ 0.45 cV in 
Fig. [21 In addition, a longer groove can result in more 
pronounced and sharper peaks; this is evident from the 
absorbance of ft. = 2.0 /J,m. In such a case, A has strong 
and non-dispersive peaks at A ~ 3.2 and 2.0 fim, respec- 
tively. The peaks at A 2.0 ^m are very close to 1 and 
are highly non-dispersive. The absorbance peak drops 
~ 10% only when the incident angle changes from 0° to 
60°, and the peak position remains close to A ~ 2 /itm. 

We now compare the TE thermal emission spectra of 
different gratings, and the results are shown in Fig. H) 
The thermal emission spectra are calculated by assuming 
Kirchhoff's law 0], and the thermal emission spectrum 
u\i{X,T) of wavevector k at temperature T is given by 



Uk(A,T) -Mb,k(A,r) X Ak(A), 



(2) 



where 'Ub,k(A, T) is the Planck spectrum of blackbody 
radiation. As our primary aim is to investigate the ef- 
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FIG. 4: Absorbance spectra of tungsten gratings with dif- 
ferent incident angles: 9 = 0° (lines) and 9 = 60° (lines with 
symbols). The normal absorbance spectrum of a 0.5 fim thick 
tungsten slab (blue dashed line) is shown as a reference. The 
apparent minor irregularity in the dispersion can be traced 
to corresponding irregularities in experimental ew- The ab- 
sorbance peaks at A ~ 2.8 fim and A ~ 3.2, 2.0 /xm correspond 
to the peaks at - 0.45 eV in Fig. [li and E - 0.39, 0.62 eV 
in Fig. Eli. 



at infrared frequencies (A > 2 ^m), which is a direct 
consequence of the high level of absorbance (Figs. [2][3]). 
Similar emission characteristics are observed when the 
temperature is raised to 2000 K (Fig. [5]). The emission 
is now dominated by photons with shorter wavelengths 
(1.0 /im < A < 1.5 /im), with the exception of A 2.0 
/im for h — 2.0 /im. In other words, the thermal emission 
of photons at A 2.0 /im is highly robust in all angles 
and a wide range of temperatures (^ 500 K). 
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FIG. 5: Thermal emission spectra of the tungsten gratings 
at 0° (lines) and 60° (lines with symbols) at 1500 K. The cor- 
responding flat slab spectrum (blue dashed line) and Planck 
spectrum (red bold line) are also shown. 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 



h = 0.50, 6 = 

^ /! = 0.50, 6 = 60 

/i = 2.00, 6 = 
<s. /i = 2.00, 6 = 60 

slab 

Planck 




2.0 3.0 
X (|im) 

FIG. 6: Thermal emission spectra at 2000 K. 



feet of the groove on the radiation spectra, we ignore 
the temperature effect on the tungsten dielectric constant 

MM- 

In Fig. [5l the thermal emission spectra of the surfaces 
corresponding to those in Fig. [4] at 1500 K are drawn 
and compared with that of a blackbody (wpianck)- In all 
cases, the emission of the grating in all frequency ranges 
is enhanced with respect to that of a slab, with the cor- 
responding optical spectra being close to that of a black- 
body. Also, the peaks of the grating's emission spectra 
are effectively red-shifted due to the strong absorption 



The results shown in Figs. suggest that gratings 
are effective in manipulating the thermal emission of pho- 
tons. By controlling groove depth, the emission peak can 
be tuned to the desired frequency. The reason for such 
emission enhancement could be attributed to the excita- 
tion of the waveguide mode, in which thermally excited 
photons at the resonant frequencies accumulate in the 
vicinity of the groove. As the excited field is strong in 
the air, but weak on the metallic surface, energy loss 
is minimal, and photons can be emitted. To show the 
presence of the waveguide modes, we consider, for Fabry- 
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Perot-like resonance, the wavelengths of the ith and jth 
resonances are related by 



According to Figs. [S]and[Sl A w 3.2 /xm and 2.0 /xm at res- 
onance, which gives ? = 1, j = 2, (2j -|- l)/(2i -f 1) = 1.66 
and A1/A2 = 1.6 « 1.66. As the waveguide mode is 
highly frequency-selective, those thermally excited pho- 
tons that are not at the resonant frequencies are re- 
absorbed by the groove and cannot be emitted. 

It should be noted that the emission peak can be 
tuned according to the operating temperature. For in- 
stance, the emission peak of blackbody radiation is given 
by Wien's displacement law, Amax = b/T, where b = 
2.897 X 10"'^ m-K is the Wien's displacement constant. 
For T — 1500 K, Amax ^ 1.93 ^m, which is close to 
the absorption peak of the h = 2.0 /Ltm grating shown in 
Fig. [3] Therefore, absorption spectra (Figs. [2]|4]) provide 
useful information on how to achieve maximum emission 
at the desired temperature. 

In this study, the thermal property of emission is mod- 
ified by enhancing the emitting power of photons at fre- 
quencies at which the emitting power would be low for 
a plain slab. Another way to manipulate thermal emis- 
sion is to suppress the emission of photons by introducing 
photonic or plasmonic band gaps to the system. The ulti- 
mate goal for thermal emission control is to acheive both 
of these on the same geometry, that is, to enhance emis- 
sion in the desired frequency range while suppressing the 
emission of un-wanted photons. Also, in practice, the in- 
trinsic thermal properties should be taken into account. 
For example, tungsten is used in this work because of its 
high melting temperature (~ 3700 K). Depending on the 
operating temperature, one may choose an emitter made 
of other metals, and the resultant emission properties 
will be different from those in this work. Nevertheless, 
the principle of emission manipulation should remain the 
same. 

The present study focuses on the effect of groove depth. 
However, the period of the grating is also an important 



parameter for controlling its optical properties. For ex- 
ample, the excitation of the SPPs shown in Figs. is 
mainly controlled by periodicity, and the excitation of 
SPPs via various geometries has been actively studied. 
In partcular, it has been shown that [s^, [36j the thermal 
radiation emitted via SPP should be highly directional 
and coherent, and will resemble the state of coherence 
of laser radiation. The tungsten gratings are thus ex- 
pected to have SPP-originated emission properties simi- 
lar to those of the SiC grating of Greffet et al. [IBl . How- 
ever, the absorbance spectrum shown in Fig. [3^ is fairly 
broadband even when SPPs are excited aX E > 0.6 eV 
(see the green dashed lines in Fig. [3^), and the emis- 
sion is merely driven by the intrinsic absorption of bulk 
tungsten. This should be due to the small periodicity 
used (flQ = 1.0 /im) in this work (compared to that of ~ 
6.25 /im in the work of Greffet et al. [3^). One should 
be able to observe SPP-driven, directional emission when 
the periodicity is increased, which will eventually lower 
the exicitation frequencies of SPP modes to frequencies 
below 0.6 eV. Nevertheless, the role of SPPs in the ther- 
mal emission of grating will be the subject of further 
study. 

Finally, the calculated absorbance spectra in this work 
can be verified in the laboratory by measuring the zcroth 
order angle-dependence reflectance. As the frequencies of 
interest are lower than the lowest SPP excitation frequen- 
cies, the higher-order diffraction effect can be neglected. 
However, the measurement of emission at high tempera- 
tures remains a challenge. To conclude, our results sug- 
gest that thermal emission can be controlled flexibly via 
metallic arrays of different groove depths, thus resulting 
in an efficient infrared radiator. It is hoped that this 
work will provide useful insights into the design of tun- 
able thermal emitting devices. 
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project no. 403308). 



[1] M. Planck, The Theory of Heat Radiation (Dover, New 
York, 1959). 

[2] R. Siegel and J. R. Howell, Thermal Radiation Heat 
Transfer (Taylor and Francis, New York, 2002). 

[3] J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Pho- 
tonic Crystals: Molding the Flow of Light (Princeton 
Universiry Press, Princeton, N. J., 1995). 

[4] K. Sakoda, Optical Properties of Photonic Crystals 
(Springer, Berlin, 2001). 

[5] H. Raether, Surface Plasmons (Springer, Berlin, 1988). 

[6] W. Y. Zhang, X. Y. Lei, Z. L. Wang, D. G. Zheng, W. Y. 
Tarn, C. T. Chan, and P. Sheng, Phys. Rev. Lett. 84, 
2853 (2000). 

[7] W. Y. Zhang, C. T. Chan, and P. Sheng, Opt. Express 



8, 203 (2001). 

[8] V. Yannopapas, N. Stefanou, and A. Modinos, Phys. Rev. 

Lett. 86, 4811 (2001). 
[9] S. Strauf, N. G. Stohz, M. T. Rakher, L. A. Coldren, 

P. M. Petroff, and D. Bouwmeester, Nat. Photon. 1, 704 

(2007). 

[10] C. Billaudeau, S. Collin, F. Pardo, N. Bardou, and J.-L. 

Pelouard, Appl. Phys. Lett. 92, 041111 (2008). 
[11] P. J. Hesketh, J. N. Zemel, and B. Gebhart, Nature 324, 

549 (1986). 

[12] C. M. Cornelius and J. P. Bowling, Phys. Rev. A 59, 
4736 (1999). 

[13] J. B. Pendry, J. Phys.: Condens. Mat. 11, 6621 (1999). 
[14] J. G. Fleming, S. Y. Lin, I. El-Kady, R. Biswas, and 



5 



K. M. Ho, Nature 417, 52 (2002). 
[15] J. T. K. Wan and C. T. Chan, Appl. Phys. Lett. 89, 
041915 (2006). 

[16] S. E. Han, A. Stein, and D. J. Norris, Phys. Rev. Lett. 

99, 053906 (2007). 
[17] J. C. W. Lee and C. T. Chan, AppL Phys. Lett. 90, 

051912 (2007). 

[18] J.-H. Lee, W. Leung, T. G. Kim, K. Constant, and K.-M. 

Ho, Opt. Express 16, 8742 (2008). 
[19] J.-H. Lee, J. C. W. Lee, W. Leung, M. Li, K. Constant, 

C. T. Chan, and K.-M. Ho, Adv. Mater. 20, 3244 (2008). 
[20] C.-M. Wang, Y.-C. Chang, M.-W. Tsai, Y.-H. Ye, C.-Y. 

Chen, Y.-W. Jiang, Y.-T. Chang, S.-C. Lee, and D. P. 

Tsai, Opt. Express 15, 14673 (2007). 
[21] K. Ikeda, H. T. Miyazaki, T. Kasaya, K. Yamamoto, 

Y. Inoue, K. Fujimura, T. Kanakugi, M. Okada, 

K. Hatade, and S. Kitagawa, Appl. Phys. Lett. 92, 

021117 (2008). 

[22] O. D. Velev and E. W. Kaler, Adv. Mater. 12, 531 (2000). 
[23] P. Jiang, J. Bertone, K. Hwang, and V. Colvin, Chem. 

Mater. 11, 2132 (1999). 
[24] G. von Frcyrnann, S. John, M. Schulz-Dobrick, E. Vekris, 

N. Tetreault, S. Wong, V. Kitaev, and G. A. Ozin, Appl. 

Phys. Lett. 84, 224 (2004). 



[25] V. Yannopapas, Phys. Rev. B 73, 113108 (2006). 

[26] K. Joulaiii, J. -P. Mulct, F. Marquicr, R. Cariniiiati, and 

J.-J. Grcffot, Surf. Sci. Rep. 57, 59 (2005). 
[27] C. J. Fu, Z. M. Zhang, and D. B. Tanner, Opt. Lett. 30, 

1873 (2005). 

[28] C. Luo, A. Narayanaswamy, G. Chen, and J. D. 

Joannopoulos, Phys. Rev. Lett. 93, 213905 (2004). 
[29] P. M. Boll, J. B. Pendry, L. M. Moreno, and A. J. Ward, 

Cornput. Phys. Commun. 85, 306 (1995). 
[30] E. RophaeU and S. Fan, Appl. Phys. Lett. 92, 211107 

(2008). 

[31] E. D. Palik, ed., Handbook of Optical Constants of Solids 

(Academic Press, Orlando, 1985). 
[32] J. A. Porto, F. J. Garcia- Vidal, and J. B. Pendry, Phys. 

Rev. Lett. 83, 2845 (1999). 
[33] G. S. Arnold, Appl. Opt. 23, 1434 (1984). 
[34] L. S. Dubrovinsky and S. K. Saxena, Phys. Chem. Miner. 

24, 547 (1997). 

[35] J.-J. Greffet, R. Carminati, K. Joulain, J.-P. Mulet, 
S. Mainguy, and Y. Chen, Nature 416, 61 (2002). 

[36] V. Yannopapas and N. V. Vitanov, Phys. Rev. Lett. 99, 
053901 (2007). 



